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Abstract. A model with a scalar type eight-fermion interaction is investigated in 
curved spacetime. The ground state of the model can be obtained by observing the 
effective potential. Applying the Riemann normal coordinate expansion, we calculate 
an effective potential of the model in a weakly curved spacetime. The result is extended 
to models with multi-fermion interactions. We numerically show the behavior of the 
effective potential and find the phase structure of the model. 
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1. Introduction 

It is expected tliat a fundamental theory with a higher symmetry may be realized at 
high energy scale. The theory is inspected in critical phenomena at early universe. A 
mechanism to break the symmetry has played a decisive role to construct a theory 
of particle physics. One of possible mechanisms is found in hadron physics. A 
composite operator constructed by quark and anti-quark develops a non-vanishing 
expectation value and the approximate chiral symmetry is broken dynamically. The 
broken symmetry is restored at high temperature, density and/or strong curvature. 
These phenomena are understood by non-perturbative dynamics in non-Abelian gauge 
theories (QCD, SM, GUT, ■ ■ ■)• 

A four-fermion interaction model is introduced to particle physics by Nambu and 
Jona-Lasinio to study low energy phenomena of the strong inter action [1]. The model 
has a similar apparent symmetry to the low energy limit of QCD. Non-perturbative 
phenomena in the model can be evaluated trough the expansion in terms of the number 
of fermion flavors. It can be observed that the chiral symmetry is dynamically broken 
for a sufficiently strong coupling according to a vacuum condensate of the composite 
field of fermions. Now some kinds of four-fermion interactions are used to analyze phase 
transition and critical behavior under extreme conditions. 

There are many works to study the curvature induced phase transition in four- 
fermion interaction models. The pioneering works have been done by Itoyama, 
Buchbinder and Kirillova in two dimensions. [21 [3j In Refs. ^ E] and [6], four-fermion 
interaction models are studied in four, three and arbitral dimensions {2 < D < A) 
respectively. The phase structure of the model is found in the maximally symmetric 
spacetime [71 [8] , Einstein space [9] and a curved spacetime with a nontrivial topology 
[To] . For a review, see for example Ref . [11] . The phase structure of the models depends 
on the coupling constants and the spacetime curvature. The broken chiral symmetry is 
restored if the spacetime curvature is positive and strong enough. However, the chiral 
symmetry is always broken down in a negative curvature spacetime even if the four- 
fermion coupling is small. This is one of the characteristic features in the four-fermion 
models. 

Since the four-fermion interactions are the lowest dimensional operator to describe 
the fermion self interaction, we often consider only the four-fermion interactions at low 
energy. However, it is not always valid to neglect higher dimensional operators in the low 
energy effective models. For example, 't Hooft introduced a determinantal interactions 
in a low energy effective model of QCD to deal with the explicit breaking of the Ua{^) 
symmetry [12]. R. Alkofer and I. Zahed considered an eight-fermion interaction to 
explain the pseudoscalar nonet mass spectrum [13] . An influence of higher derivatives 
is considered in Ref. [H]. 

In the present paper we consider higher dimensional scalar type interactions and 
study the curvature induced chiral symmetry breaking. In Sec. 2 we investigate a model 
with four- and eight-fermion interactions in curved spacetime. Using the auxiliary 
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field method and applying the Riemann normal coordinate expansion, we calculate 
the effective potential in a weakly curved spacetime. In Sec. 3 we introduce multi- 
fermion interactions and generalize the auxiliary field method. In some specific cases 
we show a correspondence with the mean field approximation. In Sec. 4 we study the 
form of the effective potential near the origin. It is shown that the chiral symmetry 
is always broken in a negative curvature spacetime. In Sec. 5 we numerically evaluate 
the effective potential and show the phase structure of the model with four- and eight- 
fermion interactions. We discuss the contribution from the higher dimensional operator 
to the chiral symmetry breaking. Finally we give some concluding remarks. 



2. Eight-fermion interaction model in curved spacetime 



In this section we introduce a model with scalar type four- and eight-fermion interactions 
and calculate the effective potential of the model in weakly curved spacetime. 

Four-fermion interaction models are often used to study the critical phenomena of 
dynamical symmetry breaking. The simplest model to induce the dynamical symmetry 
breaking is Gross-Neveu model[15|. We extend the action of the Gross- Neveu model 
with an eight-fermion interaction, 
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(1) 



where the index i shows the flavors of the fermion field ip, N is the number of flavors, 
Gi and G2 the coupling constants for the four- and the eight-fermion interactions 
respectively, g the determinant of the metric tensor, 7''(x) the Dirac matrix in curved 
spacetime and V^ip the covariant derivative of the fermion field ip. 

The action ([1]) is symmetric under the discrete chiral transformations. 



This chiral symmetry prevents the action from having mass term. 
The action is also invariant under a global flavor transformation, 



(2) 



31 



(3) 



where are generators of the flavor SU{N) symmetry. The flavor symmetry allows 
us to work in a scheme of the expansion. Below we neglect the flavor index for 
simplicity. 

For practical calculations it is more convenient to introduce auxiliary fields ai and 
CT2 and rewrite the action in the following form. 



Scr 



d^x\J —g 



(4) 
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where a is defined by 



a 



1 - 



N(T2 



The equations of motion for the auxihary fields are given by 

2Gi 



ai= - 



N 
2G 



1 — —^i^, 



(J2 



2 / 7 



(5) 

(6) 
(7) 



Substituting these equations of motion into the action (jlj), we obtain the original action 

If the non-vanishing expectation value is assigned to a, a mass term for the fermion 
filed ip is dynamically generated and the chiral symmetry is eventually broken. To study 
the phase structure we want to find the ground state of the model. In the present paper 
we assume that the spacetime curved slowly and neglect terms involving derivatives of 
the metric tensor higher than second order. We also restrict ourselves to the static and 
homogeneous spacetime. In this case we can assume that the expectation values for o"i 
and (J2 are constant. 

The energy density under the constant background ai and o"2 is given by the effective 
potential. At the large A^ limit the effective potential V is obtained by evaluating vacuum 
bubble diagrams, 



+ 



+ v{a) 



where we drop the over all factor A^ and f (a) is given by 



v{a) := — Trln(x| [iY{x)V^ 



\x) 



(8) 



(9) 



where "Tr" denotes trace with respect to flavor, spinor indices and spacetime coordinate. 
The ground state should minimize this effective potential. 

Since the four- and the eight fermion interactions are nonrenormalizable, above 
f (cr) is divergent. To obtain the finite result we regularize the divergent integral by the 
cut-off method. After performing the Fourier transformation and the Wick rotation, 
we introduce the four-momentum cut-off, A. Applying the Riemann normal coordinate 
expansion, we expand the effective potential and find 



v{(t) = - 



a^A^ + A^ In 
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(10) 



We keep only terms independent of the curvature R and terms linear in R. To determine 
the ground state we can freely normalize the effective potential. Here the effective 
potential is normalized so that V^(0, 0) = 0. 
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The necessary condition for the minimum of the effective potential is given by the 
gap equations, 

cTi dv{a) da 



and 



with 
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(13) 

(T^y + (T^ (A2 + 0-2)2 J- 

Under the ground state the expectation values for ai and a2 satisfy these gap 
equations. From the gap equations (1111) and fll2p we find the following conditions for 
non- vanishing expectation values of ai and a2. 
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0, 
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(15) 



487r2 

One of the conditions (1141) and ([TBI) should be satisfied at the minimum of the effective 
potential. In the former case the expectation values for ai and o"2 are on a hyperbolic 
curve for a positive G1G2 or on an ellipse for a negative GiG2- The latter condition is 
the same as that in the four-fermion interaction model and the eight-fermion interaction 
has nothing to do with the ground state. 

At the ground state the expectation values obey the classical equations of motion. 
Substituting Eqs.([6]) and ([7]) into the conditions (IMI) or (|T5l) . we find one of the following 
relationships, 



(i^P^f) = {i^^P)y , or {i^^py)=0. (16) 

These relationships mean that our result at the leading order of the 1/A^ expansion 
coincides with the one in the mean field approximation except for the case, (('?/''?/')^) = 0. 

We numerically evaluate the effective potential ([8]) and show the phase structure of 
the model in Sec. 4. 



Multi-fermion interaction models in curved spacetime 



6 



3. Multi-fermion interactions 



The results in the previous section can be extended more general cases. The extension 
of the model is not unique. Here we assume that vector, tensor type interactions and 
interactions with derivative operator do not develop expectation values at the ground 
state and contribution to the phase structure is negligible. We consider models with 
only scalar-type multi-fermion interactions and calculate the gap equation. 
We start from the action defined by 



S 



d^x\J —g 



2k 



k=l 



;i7) 



This action is invariant under the discrete chiral transformation ([2]) and the global flavor 
transformation ([3]). 

According to the functional integral formalism, the generating functional is given 

by 



where we set the path- integral measure, if) 
covariance. 

We consider a Gaussian integral 



-gip and ip = y^—gi', to keep the general 
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k=l 



X 



exp J d'^xV~9 



N 



1 ' 



771=1 



,(19) 



and inserts it in the right-hand side of Eq. fllSp . Thus the generating functional (1181) is 
rewritten as 

n 



c 



(20) 



k=l 



where the action Sy is given by 



S,, 



d^x^ —g 



n ^ 
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(21) 



k=l \ 1=1 

Since the number of arbitrary parameters, aki, is n{n + l)/2, we can choose a^/ 
aki{ak+i, crfc+2, cr„) to satisfy 

n ^ / ^ k \ 2 



fe=l 



1=1 
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= ficr)^^ + Yl + ' (22) 

where f{a) is a function of {k = 1, 2, ■ ■ ■ , n). Therefore the multi-fermion interaction 
terms in Eq. fl2T]) are canceled out and the action 5*^ reduces to 



d^X\J—g 



N 



k=l 



(23) 



Applying the similar method mentioned in the previous section, we obtain the effective 
potential 

" 2 

F(ai,---,a„) = 5^^ + t;(/(a)), 



fc=i 



(24) 



where v is defined in Eq.l ITOl) . 

We calculate the explicit expression of /(o") for some cases. For n = 2 it is easy to 
reproduce the result in Eq. ([5]). As another simple case, we consider the action. 



S 



d'^xyj—g 



^^7'^(a;)V^^ + 5^|^(^^) 

k=l 

In this case the function /(a) is given by the following induction series 
f[a) := (tJ,(T^^ := 



(25) 



9k 



(/c = 1, ■ ■ ■n),a^ 



0. 



(26) 



It is straightforward to show that above /(a) satisfies, 
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(27) 
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(29) 



for k = 1, ■ ■ -n. It is a sufficient condition for Eq. (!22!) . Hence, the action (p5|) is simplifies 

to 

k=i 

Since the effective potential has the same form with Eq. (l24l) . the gap equations for 
the action fl29l) are given by 

9r(cri,---,cr„) o-fc dv{f{a))df{a) 



dak 



+ 



2gk df{a) dak 



0, (A; = 1,2, 



,n] 



(30) 
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Differentiating /(a) with respect to 0"^, we get 



dfia)_^d^da[ da^^da^Jy^ -N/g, _ Naj^, 

dak da[da', da'^ da^ 

* V 9i 

for A; = 1, 2, ■ ■ ■ , n. We calculate the ratio between {k — l)-th and fc-th equations of (150]) . 

9k-i cTfc _ / df{a) \ I / dj{o) \ ^^^^ 

Substituting the Eq. (!3Tl) into Eq. (!32l) . we obtain the following configuration for ak at 
the ground state, 

1 ^^^Li 1 . ...^ 

^2 ^fc = 0- (33) 



2gkJ. 4(7^_i 1 



9fc 



9k-l 



For the action (l29l) the equations of motion are given by 



<^k = — rr\l W • 34 

^ V 9k 

At the ground state it is valid to insert Eq. (l511) into Eq. (l551) . Then we eliminate ak-i, 
CTfc, cr^ and o^'|^_^_l from both the sides in Eq. fl33l) and find following relationships between 
expectation values, 

mf) = {,p^pf, and((V^^f' )=0. (35) 

for k = 1,2, ■■■,/c' and k' G {2, ■ ■ ■ , n + 1}. The number k' depends on details of 
the effective potential and is fixed by observing the minimum of that. Therefore the 
solution of the gap equations (l30il are equivalent to the ones obtained by the mean fields 
approximation for k < k' . 



4. Chiral symmetry breaking in a negative curvature spacetime 

In a scalar type four-fermion interaction model the chiral symmetry is always broken in 
a negative curvature spacetime. Here we study the behavior of the effective potential 
near the origin and show the chiral symmetry breaking for i? < in a multi-fermion 
interaction model. 

First we consider the eight-fermion interaction model ([1]). To obtain a precise form 
of the effective potential we evaluate the curvature of the effective potential at the limit 
o"! — i> and ^ 0. Differentiating the effective potential ([8]) twice and using the 
condition (!22l) . we get 

1 - ^ 37^' (36) 



daf 2Gi V Gi J da^ 

and 

d^V 1 Na2 fldv d%\ 

+ "TT • (3 ' ) 



dal 2G2 2G2 \ada da^ 
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For a large cTj the effective potential ([H]) is almost proportional to af. Thus the potential 
has a stable minimum. If d'^V/daf or d'^V/daf is negative at the origin, the effective 
potential is unstable and the chiral symmetry has to be broken. Taking the limit ai ^ 
and (J2 — 0, we obtain, 



and 



daf 



lim 

cr^O 



1 



+ 



R 



2Gi 47r2 47r2 



In 



A2 



(38) 



(39) 



The right-hand side in Eg. (1381) is always negative for i? < 0. It implies that the auxiliary 
fields, (Ti and ct2, develops non- vanishing values at the minimum of the effective potential. 
Thus only the broken phase realizes in a negative curvature spacetime. 

Next we consider the multi-fermion interaction model f|T7|l . Differentiating both 
the sides of Eq.( !22l) in terms of (Xi, we obtain 

dficy) 1 



-aii(cr2, (73, ■ ■ ■ , (T„ 



(40) 



dai 2Gi ' 

It should be noted that that the function an does not depend on ai by definition. Using 
the condition (HOl) . we differentiate the effective potential (12^ twice in terms of cxi and 
get 



+ 



an 



dai^ 2Gi ' AGi^da^' 
To find the behavior of the effective potential near the origin, we take the limit 
(72, o"3, ■ ■ ■ , (T„ ^ 0. At the limit Eq. (!22|) we obtain, 

aii(cT2,a3,---,a„-*0) = 2Gi. (42) 

Taking the same limit in Eq.lHTl) and substituting Eq. (H2|) . we obtain Eq.( l38i) again. 
According to the same argument as the eight-fermion interaction model at least one of 
the auxiliary fields, ui, develops a non- vanishing value at the minimum of the effective 
potential. Therefore we conclude that only the broken phase realizes for the scalar type 
multi-fermion interaction models in a negative curvature spacetime. 



5. Numerical analysis of the phase structures 

Here we study the contribution from higher dimensional operators to chiral symmetry 
breaking in curved spacetime. We consider the model with four- and eight-fermion 
interactions and numerically evaluate the effective potential ([8]) with (ITOl) . We normalize 
the mass scale by the cut-off scale A. The auxiliary field a plays a role as an order 
parameter for the chiral symmetry breaking. 

We start to calculate the effective potential with fixed coupling constants Gi and 
G2- For G2 = the model contains only the four-fermion interaction. In this case the 
chiral symmetry is broken down for GiA^ > 211'^ = G^^~^^ in Mikowski spacetime, i.e. 
R=0. In a negative curvature spacetime we observe only the broken phase. |1] 
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For a model with a finite G2 we numerically evaluate the effective potential under 
the constraints f|T4|) and f|T5l) . We observe that the minimum of the effective potential 
satisfies the constraint ( fT4|) in all cases we consider here. Therefore we only show the 
results under the constraint ([HD below. 
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Figure 1. Behavior of tlie 

effective potential ® for GiA^ = Figure 2. Behavior of the 

20 and G'2A^ = 1000 effective potential ^ for GiA^ = 

20 and G2A^ = 8500. 

First we consider the model with a four-fermion coupling, Gi{> G'[f~°''), larger 
than the critical one for R = and G2 = 0. In FigsIT]and[2]we draw the behavior of the 
effective potential as a function of ^2 for GiA^ = 20, GgA*^ = 1000 and GiA^ = 20, 
G2A^ = 8500 respectively. In FigiH the chiral symmetry is broken for a negative 
curvature and is restored through the first order phase transition at i? = Rcr, as the 
curvature increases. Then the symmetric phase is realized for R > Rcr- In Fig|2]we plot 
the effective potential for a larger eight-fermion coupling, 6*2- As is seen in the figure, 
the broken chiral symmetry is restored through the first order phase transition. The 
critical curvature for G2A^ = 8500 is larger than that for G2A® = 1000. Thus the chiral 
symmetry breaking is enhanced by the eight-fermion interaction. 

Observing the minimum of the effective potential, we obtain the expectation value 
((T2) at the ground state. The expectation value satisfies the gap equations. A non- 
vanishing (cTi) is obtained for a non- vanishing and a negative ((T2) from the constraint 
( |T^ . It implies a non- vanishing (a) which corresponds to the dynamically generated 
fermion mass. Therefore the chiral symmetry is broken down for a non-vanishing and 
a negative (o"2). On the other hand the expectation value (a) is vanish for {(72) = 0. 
In this case the dynamical fermion mass is not generated and the chiral symmetry is 
guaranteed. Therefore we can inspect the chiral symmetry breaking by the expectation 
value ((T2). 

In FigsJS] and H] we draw the expectation value ((12) at GiA^ = 20 as a function 
of the spacetime curvature, R. It is clearly seen that the behavior of each lines are 
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Figure 3. Behavior of the dynam- 
ically generated fermion mass for 
GiA2 = 20 and GsA^ = 500 - 
2000. 



Figure 4. Behavior of the dynam- 
ically generated fermion mass for 
GiA^ = 20 and GsA^ = 8000 ~ 
9500. 



consistent with FigsH] and [2J The expectation value develops a non- vanishing value for 
a smaller R and disappears, as the curvature increases. We observe a mass gap in Figsl3] 
and m which reflects the nature of the first order phase transition. 

Next we evaluate the model with a four-fermion coupling, Gi{< G[^~'^^), smaller 
than the critical one for R = and G2 = 0- An interesting behavior is found for 
a smaller GiA^. At a glance we may observe a symmetric phase even for a negative 
curvature spacetime at GiA^ = 15 in FigjSl 
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Figure 5. Behavior of the 
effective potential dH) for GiA^ = 
15 and G2A8 = 8500. 



Figure 6. Behavior of the dynam- 
ically generated fermion mass for 
GiA2 = 15. 



To find the precise behavior of the effective potential we solve the gap equation. 
In Figj6] we plot the solution of the gap equation as a function of the curvature. We 



Multi-fermion interaction models in curved spacetime 



12 



observe that the expectation value, ((T2), smoothly disappears as the curvature increases. 
The transition from the broken phase to the symmetric phase is crossover. Therefore 
it is difficult to find the critical curvature of these figures. As is shown in Sec. 4, the 
expectation value ((T2) develops a non- vanishing value and the chiral symmetry is always 
broken for R < 0. 




Figure 7. Behavior of the 
effective potential (|8]) for GiA^ = 
1000 and GaA^ = 20000. 



Figure 8. Behavior of the dynam- 
icaUy generated fermion mass for 
GiA2 = 1000. 



We also consider much larger GiA^. In Fig|7] the effective potential is drawn for 
GiA^ = 1000 and G2^^ = 20000, although the coupling constant Gi seems to be too 
large as a reahstic model. It should be noted that we do not use a perturbative expansion 
in terms of the coupling constant. We observe four extrema at R/A^ = —60. As is shown 
in FiglHl we find the two kinds of the gaps for {0-2) at a negative R and a positive R. As 
increasing the curvature, R, the chiral symmetry breaking is much enhanced at the gap 
with a negative R and the broken chiral symmetry is restored through the first order 
phase transition at the gap with a positive R. The gaps appear at a large curvature 
|i?/(A^)| > 1. Since we work in the Riemann normal coordinate expansion and keep 
only terms up to linear in R, we may spoil the validity of the expansion. These types 
of the phase transition in Figs. [7] and [8] may be modified by a higher order correction 
about R. 

By observing the behavior of {o'2), we obtain the phase boundary which separates 
the chiral symmetric and the broken phases. The boundary curves are shown in FigjO] 
Below the lines the chiral symmetry is broken. The first order phase transition is shown 
by solid lines. The crossover boundary is drown by dashed lines. We find that the chiral 
symmetry is broken down for an eight-fermion coupling large enough. In the case of the 
crossover the critical value of the curvature, Rcr, is almost independent on the eight- 
fermion coupling, G2. If the phase transition is of the first order, the chiral symmetry 
breaking is enhanced for a larger eight-fermion coupling G2- 

In this section we confine our selves in the eight-fermion interaction model and 
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G2A* 

Figure 9. Phase diagram of the 
eight-fcrmion interaction model 
© for GiA2 = 15, 20 and 25. 

numerically evaluate the effective potential and solve the gap equation. It is possible to 
apply the procedure used here to the multi-fermion interaction models. Since the multi- 
fermion interaction models have more parameters, we expect that a more complex phase 
structure is realizes. Before analyze the phase structure more general models with many 
parameters, it is better to fix the physical phenomena to apply the model and reduce 
the parameters. 

6. Conclusions 

We have considered the multi-fermion interaction models as effective models to describe 
dynamical symmetry breaking at high energy scale and investigated the phase structure 
in weakly curved spacetime. Applying the Riemann normal coordinate expansion, we 
have obtained the explicit expression for the effective potential in the leading order of 
the 1/N expansion. One of the solutions of the gap equation coincides with that in the 
mean field approximation. 

To inspect the existence of a non-trivial ground state we have analytically evaluated 
the second derivative of the effective potential at the trivial solution of the gap equation, 
o"! = (72 = 0. The result is always negative for R < 0. It implies that the state with 
o"i = = is unstable. Thus only the broken phase can be realized in a negative 
curvature spacetime. We conclude that the chiral symmetry is always broken in a 
negative curvature spacetime at the large limit . It is one of the characteristic 
features of the multi-fermion interaction model. 

We have numerically evaluated details of the phase structure for the eight-fermion 
interaction model. The broken chiral symmetry was restored at a certain critical 
curvature, -Rcr(> 0). The contribution of the eight-fermion coupling enhances the chiral 
symmetry breaking. For a small four-fermion coupling, Gi < G^^~^\ we observe the 
crossover behavior, i.e. the non- vanishing expectation value, (a), smoothly disappears. 
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as the curvature R increases. For Gi > G[^~'^^ the mass gap appears for the solution of 
the gap equation and only the first order phase transition realized. We found the phase 
boundary dividing the symmetric phase and the broken phase. There is no symmetric 
phase for a negative R, as is shown analytically. 

We also apply our result to a strongly curve spacetime, though we may spoil the 
validity of the weak curvature expansion. (For the four-fermion interaction model the 
validity of the weak curvature expansion in the strongly curved spacetime is discussed 
in Ref. [6j .) A more complex phase structure is found in this case. We have observed 
two local minima for the effective potential and two kinds of mass gap for a negative 
and a large positive curvature. 

We are interested in applying our results to critical phenomena in the early 
universe. A supersymmetric extension of the four-fermion interaction model is discussed 
in Ref. BIO. A contribution from supersymmetric partners may play an important role 
at high energy scale. We will continue our work further and hope to report on these 
problems. 
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